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ABSTRACT 

We report Spitzer observations of five newly identified bow shocks in the 
massive star-forming region RCW 38. Four are visible at IRAC wavelengths, the 
fifth is only visible at 24 /im. Chandra X-ray emission indicates that winds from 
the central 05.5 binary, IRS 2, have caused an outflow to the NE and SW of the 
central subcluster. The southern lobe of hot ionised gas is detected in X-rays; 
shocked gas and heated dust from the shock-front are detected with Spitzer at 
4.5 and 24 fim. The northern outflow may have initiated the present generation 
of star formation, based on the filamentary distribution of the protostars in the 
central subcluster. Further, the bow-shock driving star, YSO 129, is photo- 
evaporating a pillar of gas and dust. No point sources are identified within this 
pillar at near- to mid-IR wavelengths. 

We also report on IRAC 3.6 & 5.8//m observations of the cluster DBS2003- 
124, NE of RCW 38, where 33 candidate YSOs are identified. One star associated 
with the cluster drives a parsec-scale jet. Two HH objects associated with the 
jet are visible at IRAC and MIPS wavelengths. The jet extends over a distance 
of ~3 pc. Assuming a velocity of 100 km/s for the jet material gives an age of 
3xl0^yr, indicating that the star (and cluster) are likely to be very young, with a 
similar or possibly younger age than RCW 38, and that star formation is ongoing 
in the extended RCW 38 region. 
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Introduction 



Young stars are not only shaped by their natal environments but shape it in return. 
They are responsible for triggering further generations of star formation, ablating the cir- 
cumstellar disks of their lower mass neighbours, clearing gas and dust in parsec-scale cavities 
surrounding them, and altering the processing of dust grains in the intercluster medium. Ob- 
servations of young clusters in the infrared (IR) and X-ray have been combined to allow the 
identificati on of the stellar popu lation in these regions from the most massive down to the 
substellar ( jPreibisch et al.l 120111). Such surveys indicate that perhaps more th an 60% of 
young stars form in clusters (jCartwright fc WhitworthI 12004 : iMegeath et al.ll2004j ). Our own 
Sun is thought to have formed in a massive cluster and have assi milated material (s uch as 
^^Fe and other short-live nucleotides) from a near-by supernova (iHester et al.l 120041 ). The 
study of the effects of high mass stars in star forming regions is thus of great importance in 
the understanding of both the large scale development of clusters and of the smaller scale 
evolution of planetary systems around lower mass Sun-like stars. Such studies indicate that 
the sub-structure within these sites of star formation is complex and this can have an impact 
on the natal environment of individual stars. 

Nearby regions of high mass star formation have recently been shown to exhibit varying 
spatial morphologies, such as the corapact spherical structure of Tr 15 (IWang et al.ll201l[ ) 



or the subclustered Tr 16 (IWolk et al. 



201 ll ). where the clusters are contigous and show no 



sign of mass segregation. RCW 38 itself shows evidence for at least four subclusters which 
are more loosely associated than those of Tr 16 and show indications of mass segregation in 
the central subcluster but not in the surrounding subclusters. By studying a number of high 
mass clusters we can better understand how the spatial distribution of young stars and their 
exposure to photo-ionising radiation can affect the fomation and evolution of young stars 
and their planets. 

One way that massive stars interact with their environment is through strong stellar 
winds. These winds can clear gas from their environs and in so doing generate shocks where 
the high velocity outflowing gas intersects wit h the ambient interstell ar /cluster medium or 
with the winds from other local massive stars ( iDraine &: McKeelll993l ). Termination shocks 
arise where the stellar wind terminates as it is balanced by the ram pressure from the 
interstellar medium (ISM) or the wind from a local massive star. If the relative velocites of 
the star and ISM/wind are large the shocked material is swept backwards around the star 
forming a 'bow shock'. The orientation of the bow shock and its 'stand-off' distance, the 
distance from the generating star to shock, are determined mainly by the the velocity of the 
driving star and its wind, and the flow of the ISM or massive star wind. 



Parsec scale collimated outflows (jets) from young protostars entrain surrounding gas 
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and dust and can disrupt their parent clouds, perhaps triggering star formation in regions 
where over densities are generated or bringing it to an end by removing the gas from which 
young stars form. Massive stars with their strong stellar winds are particularly effective 
in this regard, but Herbig AeBe and T Tauri stars are also known to drive outflows. Such 
outflows are usually identified in the visible and near-IR as Herbig-Haro (HH) objects or k nots 



of roughly symmetric emission around a YSO (IBally et al.lll995 



Ballv fc pevinelll994[) . X 



ray em ission from jets has als o been observed, e.g. in DG Tau (iGiidel et al.ll2008l ) and L1551 
IRS 5 (jSchneider et al.ll201ll ). The lifetimes of these jets are short, typically ~10^ yr. They 
are typically most active during the high accretion protostellar phase, but can c ontinue into 



the c l ass II phase. T he knots of emission may be linked to episodic accretion (IBally et al. 
1995l : lReipurthlll99lh . 



Recently, observations of cluster populations in the mid-IR with Spitzer h ave identified 



emissi on from outflows, HH-objects and bow shocks in star forming regions. iPovich et al. 



( 120081 ) report the identification using the Spitzer GLIM PSE survey da t a, of six bow shocks 
in two massive star forming regions. Ml 7 and RCW49. iTakami et al.l ( l2010l ) have modeled 
emission from six HH-objects over the IRAC bands and have found that the emission is 
mainly due to H2 emission and unlikely to have significant contributions from H I, [Fe II], 
fluorescent H2, or PAH emi ssion. The HH obj ects HH 54 and HH 7-11 have been the 
subject of a recent study by iNeufeld et al. with Spitzer IRS that mapped emission 

from rotational lines of molecular hydrogen (associated with the IRAC bandpasses) and fine 
structure transitions of ions such as Fe"*" (associated with the MIPS 24 /xm bandpass). The 
authors found that the H2 emission was consistent with nondissociative shocks with velocities 
of ~ 10 — 20 kms~^ and the fine structure emission with faster shocks with velocities of 
~ 35 — 90 kms~^ offset from the slower shocks, indicating that shocks detected at 24 /im 
may be faster dissociative shocks than those detected in the IRAC bandpasses only. 

Bow sho cks have also been consid ered as indicators for run-away 0-stars in young mas- 
sive clusters (iGvaramadze et al.l 1201 if ). The orientation of their bow shocks relative to the 
cluster, pointing away from its centre, implies the star has a velocity which is taking it out of 
the cluster. Such observations indicate that the 0-stars in question form ce ntrally and th eir 
current trajectories arise from dynamical interactions in the cluster centre ( iKroupal 120041 ). 



The RCW 38 cluster is a region of high mass star formation, lying at a distance of 
about 1.7 kpc from Earth. The region is one of the closest high mass star forming regions 
known, and contains an estimated 10^ to 10^ mem bers, including more than thirty O star 
candidates ( IWolk et al.l l2006l : IWinston et al.l 1201 ll ) . An extensive overview of the previous 
studies of RCW 38 is to be found in the Handbook of Star Formation article on the cluster 



( IWolk et al.l l2008l ) . The massive star candidates are spatially distributed throughout the 
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observed region but are concentrated in four subclusters. The central subcluster of RCW 38 
is dominated by two bright mid-IR s o urces: the 05.5 binary IRS 2 and the b ri ght emission 
ridge IRS 1 JProgel fc PerssonI Il974l : Ismith et ali IiQQqI : lOeRose et al.l 120091 ) . ISmith et al. 



( 119991 ) suggest that the central region is a cavity blown by the wind from the IRS 2 binary 
and that IRS 1 is a ridge of material swept into a shell around IRS 2. 



In our previous paper on the RCW 38 region, IWinston et al.l (120111 ) [hereafter referred to 
as Paper 1], we presented mid-IR observations taken with Spitzer IRAC &; MIPS combined 
with Chandra archival data. In this paper we will report on the identification of five bow 
shocks surrounding YSOs in the RCW 38 star forming region. We will then discuss the larger 
scale evolution of the cluster due to the winds and/or novae of massive stars in the central 
subcluster. Lastly we will discuss the embedded cluster to the NE of RCW 3 8: DB S2003- 
124, first identified in a survey for new IR clusters in 2MASS by iDutra et al.l (120031 ). and a 
newly identified protostellar jet associated with a young star in that cluster. 



Observations and Data Reduction 



The observ ations and data re duction for the Spitzer and Chandra data were presented 
in Paper 1 and IWolk et al.l ( 120061 ). A brief summary of that work is presented here. We ob- 
tained obse r vation s of the RCW 38 region with the Spitzer InfraRed Array Camera (IRAC; 
Pazio et al.l (j2004j )) in four wavelength bands: 3. 6, 4.5, 5.8, a r id 8.0 /xm and with the Multi- 
band Imaging Photometer for Spitzer (MIPS) (IRieke et al.l l2004j ). This photometry was 
supplemented by J, H and K-band photometry from the 2MASS point source catalogue 
(ISkrutskie et al.ll2006l ). The infrared observ ations were further com bined with Chandra Ad- 
vanced CCD Imaging Spectrometer (ACIS; IWeisskopf et al.l (120021 )) X-ray data. 



In Paper 1 we identified 624 young stellar objects in RCW 38, of which, 226 were 
detected in X-rays wit h Chandra. The YSOs were of the following evolutionary classes 
dlada fc Wilking||l984[ ?): 23 class I, 90 fiat spectrum, 437 class II, and 74 class III. Of the 
YSOs, 9 class I, 18 fiat spectrum, 125 class II, and 74 class III were also detected in the 
Chandra X-ray observation. The spatial distribution was studied using YSO density maps 
and minimum spanning-tree techniques and found to show evidence of three subclusters 
surro unding the central IRS2 subcluster. One of these subclusters is Obi-36 (iDutra et al. 
20031 ) associated with the reflection nebula vdBH-RN43 (Ivan de Bergh fc Herbstlll975l ) and 
has now been identified as a subcluster of the R CW 38 comp l ex. W e found that the observed 
relationship of Nh/Ak was consistent with the lVuong et al.l ( 120031 ) ratio of 1.6 x 10^^ for the 
local ISM and nearby molecular clouds. An examination of photometric variance identified 
72 variable sources, with variables in each of the subclusters and in the dispersed population. 
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Using near-IR color-magnitude criteria, 29 candidate O stars were identified. These were 
also distributed across the extended RCW 38 region, with massive stars associated with each 
of the subclusters. 

The bow shocks identified in this paper have an extended geometry such that their 
fluxes could not be obtained using the point source extraction methods applied to the stellar 
sources in Paper 1. The flux apertures were constructed manually to most accurately enclose 
the shock on the image. The area of each bow shock was estimated using an irregular 
polygonal aperture in ds9 (jjoye fc Mandell 120031 ) and the flux estimated using funtools 
( iMandel et al.ll200ll ). The same aperture was used at each IRAC band, being defined in 
the band where emission in that shock was strongest and best defined, and then replicated 
at each wavelength. The background sky emission was estimated and removed by taking 
the flux in an adjacent rectangular aperture and scaling to the area of the shock aperture. 
The 24 bow shock fluxes were estimated using funtools, similarly to the IRAC fluxes, 
using the same irregular polygonal aperture used in the IRAC bands where the shock was 
detected. 



3. Discussion 

3.1. The Bow Shocks 

In RCW 38 we identify five bow shocks in highly differing environments. The coor- 
dinates and photometry for the five driving stars are given in Table [TJ Their individual 
identifiers correspond to the Spitzer IDs listed in Table 2 of Paper 1, where their evolution- 
ary classification and variability over our two epochs of IRAC observations were determined. 
The coordinates of the apexes of the bow shocks and their estimated fluxes are listed in 
Table [21 Figure [T] presents an IRAC false-color image of the RCW 38 cluster with the loca- 
tions of the five identified shocks highlighted. The bow shocks are clearly visible in shocked 
hydrogen emission, seen predominantly at 4.5/xm (in green). PAH emission from heated dust 
is visible at 8.0/im. The region of bright yellow emission to the west of centre in the image 
is the central subcluster of RCW 38 where IRS 1 & 2 are located. Figure [2] provides IRAC 
and MIPS false-color enlarged images of each of the five identified shocks and their driving 
stars. The upper two plots (YSO 293 & YSO 49) show the bow shocks clearly visible in 
shocked hydrogen at 4.5/im. The shocks in the central plots (YSO 1 29 fc YSO 803) emi t 



predominantly at 24/zm, indicating perhaps faster dissociative shocks (INeufeld et al.ll2006l ). 
The lower plot shows YSO 581, which is seen as a faint emission feature against IRS 1, the 
bright X-shaped emission feature to the NE of the shock in the centre of the image. IRS 2 
is the bright blue source to the E of IRS 1. Details on the individual driving stars and bow 
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shocks can be found in the appendix. 



The near- and mid-IR fluxes were used to construct spectral energy distributions (SEDs) 
for the five driving stars and bow shocks, with the plots of the SEDs given in Fig. Ei The star 
symbols show the stellar fluxes while the squares indicate the fluxes of the shocks. The shock 
fluxes of YSO 803 and YSO 129 are very faint shortwards of 24/zm and so the IRAC fluxes 
are not greatly in excess of the background and are more uncertain. Both the stellar and 
shock fluxes from YSO 581 are contaminated by the superposition of the shock and star due 
to the line of sight. The objects YSO 581, YSO 293, and YSO 49 were classified as CO/I and 
II sources as their stellar SEDs exhibit excess emission above that of a stellar photosphere. 
The SEDs of the bow shocks are consistent with emission from a black body. These fi nding s 
are con sistent with those of shocks in RCW49 and M17 presented by iPovich et al.l (120081 ) . 
and by iBally et al.l (120001 ) in the Orion Nebular Cluster. The ONC is relatively nearby at a 
distance of 145 pc, M17 lies at a sim ilar distance (1.5 kp c) to RCW 38, while the distance to 
RCW 49 is given as 4.2 or 6 kpc in lPovich et al.l (120081 ). The five RCW 38 shocks increase 
the number of known mid-IR shocks by ~40% and more than double the known sample at 
1-2 kpc. 



3.1.1. Dynamical Implications 



Consideration of the stand-off distance of the bow shock (the distance from the star to 
the apex of the shock) and the mass loss rate of the driving star can l ead to an est i matio n 
of the relative velocity of the s tellar winds or ISM in t he reg ion. From lPovich et al.l (120081 ). 
following the normalisation of Ivan Buren fc McCrayl (119881 ). we have an equation relating 
the velocity (vo) and hydrogen particle density (no,3 = lO^cm"^) of the ambient medium to 
the mass loss rate (M^^_6 = MqHt'^) and relative velocity {v^fi = 10®cms~^) of the 
driving star and the bow shock stand-off distance (ci^) 



1/2 



vou,^ = 1.5(— )-i(M^,_6t;^,8)'/' 
pc 



The stand-off distance must incorporate a cos{i) term since we cannot ascertain the angle of 
the shock with respect to the plane of the sky. However, cos{i) has a mean value of 0.9 for 
angles between 0" and 45° so this will affect the resulting velocity measurement at the 10% 



level. The density of the ISM is assumed to be approximately 10" 



cm 



and hence equal to 



unity in the equation. The mass loss rates and stellar wind velocities from young massive 



stars are highly variable and strongly dependent on stellar spectral type (iFuUerton et al. 



20061 ) . This uncertainty is further complicated by the estimation of the spectral type, which 
is calculated from the near-IR colors and not spectroscopy. In this case a value of <0.2M^ _e 
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is applied to all stars with spectral type later than 09 applicab l e to a ll five objects in this 



case. This value is based on estimates taken from IVink et al.l ( 120011 ) and iFuUerton et al. 



( 2006). A thorough description of the assumptions made in this calculation is presented in 



Povich et al.l tOQ& \ 



The offset in pixels from the coordinates of the star to the apex of the bow shock was 
measured on the 8 fim mosaic and converted to distance in parsecs assuming a distance 
to ROW 38 of 1.7 kpc. The spectral types of the driving stars were estimated from their 
positions on the K v. J—K color-magnitude diagram. The spectral types, stand-off distances 
and estimates of the relative velocity for the shocks are given in Tabled except YSO 581, the 
geometry of which precludes the measurement of the stand-off distance. With the caveat that 

1/2 

the values of VouJ^ are highly uncertain we can say that in the cases of the four shocks with 
measurements the velocities span a small range, from 4-12 kms~^. The velocities calculated 
for sources 129 and 803, those only visible at longer wavelengths, are somewhat lower than 
those of the sources with shocks clearly visible in the IRAC bands. This would not agree 
with their emission arising from fast dissociative shocks and is perhaps an indication that 
these more outlying sources are located in a lower density medium. It may also be linked to 
evolutionary classification: both are class III sources. It may also be related to their being 
shocked by the ISM and not due to interaction with wind from high mass stars. 

The orientations of all five bow shocks in RCW 38 appear to be aligned with 0-star 
candidates identified in their vicinity. None of the five shocks are oriented such that their 
driving stars would appear to be following a trajectory away from IRS 2. This allows us to 
rule out the possibility that the stars driving the bow shocks are run away 0-stars that are 
being ejected from the central subcluster. YSOs 293 and 49 are aligned to ocl5-vl3 (see 
Table 11 in Paper 1), while the apex of the central bow shock is clearly pointing towards 
the central 05.5 binary IRS2. This orientation indicates that the shock is ~ 60deg SW of 
the binary and also perhaps ~ 15 — 45 deg behind the O stars. This would imply that the O 
star binary is not clearing out a cylindrical structure as this is a projection effect and these 
walls lie behind the massive stars. Their winds are likely compressing and illuminating the 
posterior wall of the H2 blister bubble. In the case of YSO 129 and YSO 803 the source of 
the interaction is less clear; it may be the candidate 0-star NE20, associated with the new 
subcluster, c.f. Sec. 13.31 However it is more likely that they trace the direction of the flow 
of the ISM material in their vicinity, implying that this material is moving in a roughly E 
to W direction across the face of the cluster. 
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3.2. Outflows & Winds: Influence of the massive stars 



Discussion of ttie central O stars of IRS 2 leads to an examination of tlie overal structure 
of RCW 38 and the effect of the massive stars and their winds on the evolution of the cluster. 
In planetary nebulae, complex X-ray emission structures have been observed which show both 
evidence of axial symmetry and an anticorrelation between the X-ray diffuse emission and 
the infrared extinction, indicating t hat the X-ray outflows can in part be shaped by the 
material into which they expand (cf: iKastner et al.ll2002l : IChu et al.l 120041 ) . Similar features 
have also been observed in the centre of RCW 38. The two 05.5 stars of IRS 2 have an 
estimat ed combined luminosity loq jL) o f 5.71 Lp, and hydrog en ionising photon flux log^Qn) 
of 49.4 JPeRose et al.lbood : ISmithlbood : iMartins et aL||2005 [) and are co nsidered likely to be 
the main drivers of the outflows from the central subcluster. IWolk et al.l ( 120021 ) modelled the 
diffuse X-ray emission and found it was best flt by synchotron emission. We infer from the 
orientation of the central bow shock of YSO 581, that the central 0-stars of IRS 2 are not 
hemmed in by the dense dust and gas structures of IRS 1 and surrounding rim of material. 
Rather, they may be some distance in front of the face of the molecular cloud, perhaps 
~0.3 pc, (the radius of a sphere centred on IRS 2 and deflned by the rim of emission at 8 /im 
at the assumed distance to RCW 38 of l.Tkpc). The outflows and stellar winds from these 
stars have been been funnelled outwards in three outflows along a NE-SW axis and to the 
NW, evacuating the region and heating the gas and dust along the path as they are swept 
up in the moving shock front. 

Further, the orientation of the outflow is possibly tilted with respect to the plane of the 
sky, with the SW lobe projecting out of the plane by a few degrees, while the NE lobe is 
projected into the molecular cloud in front of which the cluster sits. The southern lobe has 
met with less dense material at the border of the molecular cloud and ISM. This material 
is being shocked where the hot ionised gas of the expanding shell collides with it, as can be 
seen in Fig. [51 The hot gas can be observed in the X-ray contours shown in Fig. [5l fllling 
the expanding cavity cleared by the outflow. The leading shock front is visible at 4.5 i^im 
as shocked Hydrogen emission and also via the heated dust visible at 24 /im, cf. Fig. H] & 
|5l The expanding flow of hot ionised gas to the NW visible in the X-ray data has caused 
shocked emission in the IR as it collides with and flows around the dust and gas. The third 
outflow, expanding to the NE, appears to be moving into a region of higher dust density 
and appears to be partially obscured by a ridge of dense dust. These outflows may have 
triggered the latest generation of star formation in the central subcluster, as can be traced 
by the fllament of protostars, illustrated in Fig. |5l The YSOs shown in Fig. |5l identifled in 
Paper 1, are the class O/I (circles) and flat spectrum sources (triangles). While the class 
O/I protostars are located throughout the region, they are notably concentrated in discreet 
groups in the three outer sub clusters and within the central IRS 2 subcluster where they 
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form an filament or arc-shaped distribution. The protostars i n this filame nt are also spatially 



coincident with a ridge of dust emission observed at 1.2mm (j Vigil! |2004[ ) 



Outside of the central subcluster, the molecular cloud is also being shaped by other 
massive stars. A dense pillar of dust is being carved away by the photo-evaporating winds of 
YSO 129 (one of the bow shock driving stars) and possibly the 0-star candidate ovl5-vl3, 
c.f. Fig. [T] and for an expanded view Fig. O These two stars lie at projected distances from 
the top of the pillar of 1.37 and 4.94 arcmin (0.68 and 2.44 pc) for YSO 129 and ovl5-vl3, 
respectively. The pillar has a finger of dense material extending from it in the SW direction 
and two illuminated faces to SW and W. The faces and finger of the pillar show very bright 
emission at 4.5, 8.0 and 24/zm implying both heating of the dust and shocked gas emission. 
No point sources are detected in the pillar. This is perhaps because the density is currently 
too low to promote star formation or because any protostars within are as yet too young to 
be detected in the near- to mid-IR and would require observations at longer wavelengths, 
such as are available with Herschel or with ground-based submm instruments. 

Six of the candidate O stars are located on the periphery of the central subcluster. If 
these stars are early type members, then the surrounding low to intermediate mass YSOs 
that have some shielding from the IRS 2 binary would be exposed to higher fiuxes of ionising 
photons than might otherwise have been expected. This may have an impact on the X- 
ray calculated YSO disk fraction (the ratio of X-ray detected class II to class III) which 
was found to be constant with radial distance beyond ~0.5pc from the central IRS 2 binary. 
Within a radial distance of IRS 2 of 0.5pc the ratio is 40 ±28%; this rises to 67±5% beyond a 
radius of 0.5pc and is approximately constant from 0.5-3pc from the centre of the subcluster. 
While the uncertainties do overlap (due to the small sample in the central region), there is 
an indication that the disk fraction is decreased in the vicinity of IRS 2. The three other 
subclusters each contain a candidate massive star, and while not as early as IRS 2 they are 
still capable of affecting the circumstellar disks of their subcluster members. The majority of 
the candidates are strewn throughout the distributed population and thus it may be difficult 
to identify any significant fraction of the population of the extended RCW 38 region that is 
isolated from the effects of a massive star. 



3.3. DBS2003-124 & Protostellar Jet 



The embedded stellar cluster DBS2003-124 was identified by iDutra et al.l (120031 ) in a 
survey of the 2MASS catalogue designed to identify new IR clusters. Located at coordinates 
OQhvQQm^l^^s^ -47'^26"'02.3^ 16.9' NE of IRS 2, it lies close to edge of the IRAC image and 
was only observed in the 3.6 & 5.8/im bands due to the offset in the fields of view of the 
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IRAC channels. The cluster was also observed in the MIPS 24/im image, c.f. Fig. [61 Visually, 
the region appears to be contiguous with RCW 38, the cloud material and YSOs extending 
without interruption between this cluster and the subclusters associated with RCW 38 itself; 
therefore we assume the distance to also be 1.7kpc. At this distance it lies 8.4pc in projection 
from IRS 2. 

The cluster is surrounded by a ridge of bright emission at 5.8/im, within which are 
located tens of point sources, including one very bright 24^m detection in the centre of the 
grouping and another 24^m source on the NE of the central group. A string of five 24/im 
point sources is observed, two within the ridge of dust emission and a tail of three more 
sources extending to the SW, one of which appears to be the origin of the protostellar jet. 
The projected separations of these five objects is semi-regular from the northernmost source 
to the south: 0.26pc, 0.43pc, 0.27pc, 0.42pc. The southern most of the five 24/im sources, 
NE20, is a candidate high mass star, and it is located 1.16 pc from the subcluster centre. 
This would make DBS2003-124 the only one of the five identified subclusters without an O- 
star at its centre. Unlike Obj 36, the NW subcluster identified in Paper 1, there is no bright 
mid-IR emission nebula associated with this small cluster of stars. The region also exhibits 
faint extended emission at 24yum and is the only subcluster where 24/im point sources are 
identified (though both the central subcluster and Obj 36 are saturated in the 24/xm mosaic). 

Two methods were used to identify candidate members of the cluster. In the first 
method, those sources within the ridge of dust surrounding the main 24/xm source and the 
24/im sources themselves were considered to be members of the cluster and their coordinates 
and photometry are included in Table HI This method identified 14 possible members from 
the 3.6 and 5.8/tm photometry, with the five at 24/im. The second method involved searching 
for sources with excess emission on the [3.6] v. [3.6 - 5.8] color-magnitude diagram in the 
region surrounding the cluster and 24/im filament. With this method, we identify 17 objects, 
including 11 new possible cluster members. The coordinates and photometry of the sources 
identified using the cmd are presented in Table |5l This gives an estimated membership of 33 
objects in DBS2003-124. High resolution imaging of the cluster in the near-IR would allow 
for a better estimate of the total population of this young cluster. 

3.3.1. The Protostellar Jet 

Figure |6] also reveals two aligned shells of MIR emission, which appear to be shocked 
emission from two infrared candidate HH objects in a westward moving protostellar jet. 
Tracing the outflow backwards leads to the young star NE21, one of the 24/im sources in the 
filament extending from DBS2003-124. This star also may have an extended emission feature 
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at 5.8/im, which might trace the beginnings of the protostellar jet and indicating that the jet 
remains active, though this feature may be also be a data processing artifact. The coordi- 
nates of the apexes of the two shells are: 09'''•00'"13.94^-47'^27™39.05" and 08'^'^59™58.93^ 
— 47'^27™'54.28''. These are found at distances from NE21 of 1.85pc and 3.16pc, respectively, 
assuming a distance of 1.7kpc. By considering a range of velocities for the material in the jet, 
we can estimate the age of the jet from the two identified shells. In Table E l we hst the ages 
of th e jet assuming velocities of 10, 100 and 500 kms~^ (IBally et al.lll995l : iBally fc Devine 
19941 ). We estimate the age of the jet to be ~ 6 x 10'^ to 3 x 10^ years. This would imply a 



very young age for this star and the subcluster as a whole. It is possible that this subcluster 
is younger than RCW 38 itself and that star formation is progressing through the molecular 
cloud. 



4. Conclusion 

We have examined the RCW 38 cluster using mid-IR observations taken with the Spitzer 
Space Telescope's IRAC and MIPS instruments and discovered five bow shocks in the region. 
We have combined these data with 2MASS photometry and Chandra observations to clas- 
sify the evolutionary stage of the driving stars and to estimate their spectral type. We have 
examined the large scale outflow from the central subcluster to study its effects on the evolu- 
tion of the cluster as a whole, in particular its star formation history. We have also identified 
DBS2003-124 as a embedded mid-IR subcluster most likely associated with RCW 38 and 
two candidate HH-objects from a jet originating from a YSO within this subcluster. The 
main results from our study are: 

• Five bowshocks were detected by their mid-IR emission in the IRAC and MIPS band- 
passes. The shocks are associated with driving stars with estimated spectral types 
from 09-B3, based on their near-IR magnitudes and colors. Four of the driving stars 
are previously identified YSOs, with one class I, two class II and one class III objects. 
The fifth driving star was classified as being class III. Three of the bow shocks are 
clearly visible in the IRAC images, two are only clearly detected at 24 /im. 

• The relative velocities of the shocks were estimated from the stand-off distance of the 
shocks to their associated stars and their estimated mass loss rates and stellar wind 

1/2 

velocities. Though with large uncertainties, the values of VouJ^ were all found to be 
on the order of 10 km/s. Those shocks likely interacting with the ISM, YSO 129 & 
803, and only detected in the longer wavelengths, had slightly lower values. 

• The central O star binary, IRS2, is driving a powerful stellar wind that is shaping the 



- 12 - 



cluster environment. The outflow of hot ionised gas, traced by diffuse X-ray emission, 
is carving a cavity out of the molecular cloud, opening the region to the SW and 
NW where hot ionised gas flows outwards. Outflowing material to the north-east has 
possibly triggered a new generation of star formation in the region. 

• A dense pillar of gas and dust has been observed that is being photo-evaporated by the 
winds from the bow-shock driving star YSO 129 and the 0-star candidate ovl5vl3. 
No point sources are identified in the pillar at near- to mid-IR wavelengths. 

• We report on the identification in IRAC bands 3.6 & 5.8/xm of DBS2003-124, an 
IR cluster to the NE of ROW 38, as a subcluster of the RCW 38 complex. Using 
association with the ridge of bright 5.8/xm emmision and color- magnitude criteria, we 
have identified 33 sources as likely members of this cluster. 

• Two infrared candidate HH-objects from a protostellar jet, associated with NE21 in 
DBS2003-124, have been identified. Two shocked emission lobes are visible to the W 
of NE21, at distances of 1.86 and 3.16 pc from NE21. Estimating the velocity of the 
jet material gives an age of 3 x 10^ to 3 x 10^ years for the jet. 

This work is based on observations made with the Spitzer Space Telescope (PID 20127), 
which is operated by the Jet Propulsion Laboratory, California Institute of Technology under 
NASA contract 1407. Support for this work was provided by NASA through contract 1256790 
issued by JPL/Caltech. Support for the IRAC instrument was provided by NASA through 
contract 960541 issued by JPL. This publication makes use of data products from the Two 
Micron All Sky Survey, which is a joint project of the University of Massachusetts and 
the Infrared Processing and Analysis Center/California Institute of Technology, funded by 
the National Aeronautics and Space Administration and the National Science Foundation. 
This research has made use of the NASA/IPAC Infrared Science Archive, which is operated 
by the Jet Propulsion Laboratory, California Institute of Technology, under contract with 
the National Aeronautics and Space Administration. This research made use of Montage, 
funded by the National Aeronautics and Space Administration's Earth Science Technology 
Office, Computation Technologies Project, under Cooperative Agreement Number NCC5- 
626 between NASA and the California Institute of Technology. Montage is maintained by 
the NASA/IPAC Infrared Science Archive. TLB acknowledges support from NASA through 
a grant for HST program 11123 from the Space Telescope Science Institute, which is operated 
by the Association of Universities for Research in Astronomy, Incorporated, under NASA 
contract NAS5-26555. 
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5. Appendix: Details of Individual Bow Shock sources. 

YSO 581 The bow shock associated with YSO 581 is located in the heart of the 
RCW 38 cluster and likely arises from interaction with the strong winds emanating from the 
05.5 binary 1RS2. The probable driving star (YSO 581) is classified as being of evolutionary 
class 0/1. It was not detected as an X-ray source in the Chandra observation, and is not a 
candidate IR variable source. The bow shock is observed in all four IRAC bands but due 
to the high background and source density in the central cluster it was difficult to extract 
accurate values of the fiux. It does not have a MIPS detection due to the saturation of 
the the 24/im mosaic in the central subcluster. From the image, the bow shock appears 
to lie partially in front of the driving star, which would lead to the measured fiuxes of 
both objects being contaminated; thus it is possible that the identification of this object 
as having an envelope is inaccurate and that this 'envelope' is in fact emission from the 
shock. Likewise the fiuxes for the shock are contaminated by the stellar fiux, though the 
aperture used to measure the shock fiux excluded a circular aperture covering the star itself. 
Spectrosc opic data of the st ar showing absorption features would help to clarify if this is 



the case. ISmith et al.l (119991 ) have identified YSO 581 as source 'D' in their near-IR maps, 
where they mention the possibility of a cavity and ridge of emission to the west of IRS 1 
likely corresponding to emission from the bow shock. 

YSO 293 The bow shock associated with YSO 293 was found to the east of IRS 2 in an 
area not immediately affected by the winds from the central massive stars. The young star 
YSO 293 was classified as a class II member of the cluster. It is not a variable candidate, but 
was detected as an X-ray source in the Chandra observation. The source was detected with 
10 counts in the ACIS data, with Nh = 1.91 ±0.81 x 10^^ cm'^ and kT = 1.14 ±0.59 keV. 
The star is detected in the 2MASS bands, and in the three shortest wavelength IRAC bands, 
but not as a point source at 8.0 /im. The bow shock was detected at all four IRAC bands 
and at 24 fim. The bow shock and star are clearly separated on the sky in the shorter IRAC 
bands, though there is possibly some contamination at 8.0 yum and likely a larger degree of 
contamination at 24 /im. The shock likely arises at the interface of the stellar wind with 
that of the 0-star candidate ovl5-vl3 (c.f. Paper 1 Table 12). This star lies 0.88' to the 
east, or 0.45pc at the assumed distance to RCW 38 of 1.7kpc. The star is located in a region 
of elevated emission from hot ionised gas in the Chandra data where two candidate 0-stars 
(ovl5-vl3 & y244) are also found, as shown in Fig. |H 

YSO 49 The next bow shock is associated with YSO 49 and displays similar properties 
to YSO 293. This star is also to the east of IRS 2 and to the south west of YSO 293. The 
driving star, YSO 49, was classified as being class II. This classification is not affected by 
the emission from the shock. It was not detected in X-rays, though it does lie within the 
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field of view of tlie observation. This star did not show any indication of variabihty over the 
two epochs of IRAC observations. The star is detected in all 2MASS and IRAC bandpasses, 
where its stellar SED indicates excess emission in the IR. The shock is detected in the four 
IRAC bands and also in MIPS 24 fim. This star appears to be deeply ensconced in the 
surrounding molecular cloud. This shock also likely arises at the interface of the stellar 
winds from ovl5-vl3 and YSO 49. YSO 49 lies at a distance of 3.4' or 1.68pc from ovl5-vl3. 

YSO 129 The fourth bow shock is located further east of IRS 2 than YSOs 293 and 49. 
The driving star is YSO 129, which has an evolutionary classification of class III, indicating 
that it does not have an optically thick circumstellar disk, consistent with its SED in Fig. |3l 
This object was identified as a YSO from its detection in the Chandra observation. The 
source was detected with 143 counts in the ACIS data, with Nh = 0.23 ± 0.08 x 10^^ cm"^ 
and kT = 1.27 ±0.16 keV. It is not a variable candidate. The star is detected at all 2MASS 
and IRAC wavelengths, but not at 24 /im. The bow shock is visible only at 8.0 /xm and 
24 fim - indicating that it is perhaps more evolved than the others or is located in a less 
dense medium where the shocked dust is not heated t o as high a tempera ture. It may also 



indicate fine structure emission from a faster shock (INeufeld et al.ll2006l ). Measurements 
of the fluxes at 3.6, 4.5 and 5.8 fim were also obtained but are comparatively very faint. 
This star's bow shock is possibly interacting with NE20, a new candidate 0-star detected 
in the newly identified NE subcluster (c.f. Sec. 13.31) . However, while the projected 2-D 
orientation of the shock is consistent, the large distance, 3.96pc, and the topography of the 
dust structures surrounding YSO 129 would imply that the shock may be due to interactions 
with the local ISM. The star is located in a region of diffuse X-ray emission, which appears 
to show extended jet-like emission at the star, see Fig. [51 YSO 129 also illuminates a dense 
pillar of dust, which it appears to be photo-evaporating, c.f. Sec. 



YSO 803 The last of the five bow shocks is visible only at 24 and is located further 
to the south east of the central RCW 38 cluster than the other shocks, in a region at the edge 
of the illuminated cloud. The driving star was not identified as a YSO in Paper 1; its SED 
does not show any indication of excess emission in the mid-IR bands indicating that it could 
be a class III diskless member. It is outside the field of view of the Chandra observation 
hence its X-ray flux could not be measured. The star YSO 803 is detected in 2MASS and 
the three shorter IRAC bands, the emission at 8.0 /im was extended. Measurements of the 
flux of the bow shock were obtained at 5.8, 8.0 and 24 /zm, though the flux at 5.8 fim is very 
faint. The detection only at longer wavelengths may suggest that this shoc k has a higher 



veloc ity than the other four, as evidenced by fine structure emission at 24 //m (INeufeld et al. 



20061 ) . This bow shock may also be interacting with the winds from the candidate 0-star 
NE20, though the distance of 16' or 8pc means that even a small difference to the orientation 
of the bow shock will miss the star, we therefore conclude that this star is also interacting 
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with the ISM or ambient cluster medium. 
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Table 1. Coordinates and photometry of the stellar sources driving bow shocks. 





RA 


Dec 


J 


H 


K 


3.6/.ir?i 


4.5/xm, 


5.8^?Ti 


8.0/^m 


ID 


J2000 


J2000 


[mag.] 


[mag.] 


[mag.] 


[mag.] 


[mag.] 


[mag.] 


[mag.] 


293 8 


59:26.70 


-47:33:10.14 


11.696±0.024 


10.967±0.027 


10.503±0.030 


10.082±0.061 


9.716±0.101 


9.168±0.210 




49 8 


59:18.44 


-47:35:51.89 


11.104±0.023 


10.438±0.026 


10.016±0.025 


9.781±0.024 


9.555±0.035 


9.125±0.083 


7.944±0.132 


129 8 


59:56.06 


-47:33:04.4 


8.587±0.021 


8.254±0.038 


8.073±0.021 


7.999±0.009 


8.169±0.012 


7.978±0.014 


8.029±0.031 


581 8 


59:02.94 


-47:30:54.1 


11.002±0.021 


10.130±0.021 


9.607±0.019 


7.966±0.038 


6.772±0.032 


5.544±0.072 




803 8 


59:41.98 


-47:41:22.8 


9.035±0.025 


8.434±0.042 


8.186±0.030 


8.879±0.015 


8.175±0.012 


7.970±0.014 


8.027±0.021 



'^See lWinston et all ll201lh {Paper 1) for further details. 
No photometry was available for the sources at 24/im. 



Table 2. Coordinates and photometry for the bow shocks. 



Sp. 






3.6/im 


4.5/im 


5.8/im 


8.0/im 


24^m, 


ID 


J2000 


J2000 


[mjy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


293.S 


8:59:27.31 


-47:33:09.0 


16.3±1.8 


38.2±2.1 


28.9±3.8 


43.1±6.2 


<3598 


49.S 


8:59:18.62 


-47:35:45.8 


21.6±1.3 


38.9±1.6 


36.3±2.7 


68.6±4.8 


<4013 


129.S 


08:59:57.36 


-47:32:45.7 


<0.8 


<1.4 


5.4±5.0 


62.7±25.4 


1251±14 


581.S 


08:59:02.94 


-47:30:54.1 


<366 


<869 


<1127 


<6950 




803.S 


08:59:43.07 


-47:41:28.6 






6.3±2.6 


45.8±3.0 


813±16 



Table 3. Stand-off distances and estimated stellar wind properties. 



Sp. 
ID 


Est. Spectral Type 


dwcosi 
pc 


1/2/ -N-l a 

voriQ 3 (cosi) 


293 


B2-B3 


0.079 


<8.41 


49 


B1.5-B2 


0.055 


<12.19 


129 


BO 


0.177 


<3.79 


581 


B0.5-B1 






803 


BO 


0.125 


<5.37 



^Calculated for an assumed mass loss rate of <0.2M™^_6 
for each star. 
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Table 4. Coordinates and photometry for the Candidate Members of NE subcluster. 



Sn 


RA 


Dec 


J 


H 




3 .a^TTl 


5 8m 
urn 


ID 


T2ono 


T2000 


fin a c 1 








[mag. J 


1 


9:00:42.92 


-47:25:52.30 




15.22±0.11 


13.98±0.08 


12.48±0.09 


11.24±0.11 


2 


9:00:42.04 


-47:26:22.32 


16.06±0.10 


14.55±0.05 


13.84±0.05 


13.16±0.11 


12.91±0.39 


3 


9:00:43.04 


-47:26:07.38 




16.14±0.19 


14.81±0.14 


13.01±0.09 


11.85±0.18 


4 


9:00:41.42 


-47:26:05.00 


16.38±0.13 


13.59±0.06 


11.75±0.04 


9.78±0.04 




5 


9:00:39.22 


-47:26:34.58 




14.65±0.06 


12.78±0.03 


11.54±0.05 


11.02±0.12 


6 


9:00:43.09 


-47:25:39.47 


13.82±0.03 


12.12±0.03 


11.00±0.02 


9.37±0.02 


7.79±0.01 


7 


9:00:42.09 


-47:25:53.74 






13.84±0.19 


11.99±0.06 


10.55±0.08 


8 


9:00:41.24 


-47:25:54.26 






13.11±0.15 


11.04±0.12 


9.08±0.04 


9 


9:00:42.01 


-47:26:00.25 






11.43±0.05 


9.98±0.05 


8.72±0.05 


10 


9:00:39.62 


-47:26:20.94 








14.82±0.39 




11 


9:00:40.88 


-47:26:14.79 








12.00±0.13 




12 


9:00:40.86 


-47:25:59.36 








10.29±0.07 




13 


9:00:41.43 


-47:25:45.97 








13.29±0.15 




14 


9:00:39.29 


-47:25:42.83 








14.76±0.26 




15 


9:00:41.94 


-47:25:33.31 








14.19±0.21 




16 


9:00:.39.90 


-47:25:35.15 








14.90±0.38 




17 


9:00:39.31 


-47:25:36.05 








14.85±0.26 




18 


9:00:40.77 


-47:26:14.78 










10.47±0.12 


19 


9:00:41.42 


-47:26:06.34 










8.00±0.04 


20 


9:00:33.80 


-47:27:58.17 


8.40±0.02 


6.60±0.03 


5.76±0.02 




6.28±0.01 


21 


9:00:36.96 


-47:27:19.10 




15.77±0.15 


12.99±0.03 


8.93±0.01 


6.44±0.01 


22 


9:00:38.62 


-47:26:48.98 




15.64±0.13 


13.51±0.07 


11.04±0.04 


8.91±0.03 



Table 5. Coordinates and photometry for the Candidate Members of NE subcluster from 



3.6/Ltm V 3.6//m-5.8/xm Color-Color diagram. 



Sp. 
ID 




RA 
J2000 


Dec 
J2000 




J 

[mag.] 




H 
mag. 






inag. 




3.6/im 
[mag.] 


5.8^m 
[mag.] 


22 


9 


00:38.62 


-47:26:48 


98 




15 


64±0 


13 


13 


51±0 


07 


11.04±0.04 


8.91±0.03 


23 


9 


00:26.47 


-47:26:27 


64 


16.05±0.08 


14 


52±0 


05 


13 


93±0 


05 


13.46±0.12 


12.74±0.23 


24 


9 


00:26.76 


-47:26:57 


52 
















11.91d=0.06 


11.05±0.06 


25 


9 


00:29.65 


-47:27:31 


30 


13.89±0.03 


13 


19±0 


03 


12 


90±0 


04 


12.79±0.08 


12.19±0.16 


1 


9 


00:42.92 


-47:25:52 


30 




15 


22±0 


11 


13 


98±0 


08 


12.48it0.09 


11.24±0.11 


26 


9 


00:44.90 


-47:25:20 


41 










15 


08±0 


15 


12.56±0.08 


11.08±0.07 


3 


9 


00:43.04 


-47:26:07 


38 




16 


14±0 


19 


14 


SldzO 


14 


13.01±0.09 


11.85±0.18 


27 


9 


00:34.01 


-47:25:53 


45 




15 


07±0 


08 


14 


28±0 


08 


13.12±0.10 


12.20±0.18 


28 


9 


00:36.25 


-47:25:52 


88 


16.11±0.10 


14 


46±0 


05 


13 


75±0 


06 


12.57±0.07 


10.80±0.06 


29 


9 


00:38.98 


-47:26:19 


35 










14 


82±0 


13 


11.90±0.06 


10.26±0.06 


30 


9 


00:43.83 


-47:24:04 


82 




15 


87±0 


16 


14 


11±0 


06 


12.47±0.07 


11.30±0.10 


31 


9 


00:35.82 


-47:24:46 


60 


15.37±0.06 


13 


82±0 


04 


12 


98±0 


03 


12.04±0.06 


11.14±0.08 


6 


9 


00:43.09 


-47:25:39 


47 


13.82±0.03 


12 


12±0 


03 


11 


00±0 


02 


9.37±0.02 


7.79±0.01 


32 


9 


00:43.79 


-47:24:38 


94 


15.22±0.05 


12 


99±0 


03 


11 


77±0 


03 


10.43±0.03 


9.32±0.03 


7 


9 


00:42.09 


-47:25:53 


74 










13 


84±0 
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Table 6. Jet velocities and ages for the two shells. 



velocity 


velocity 


1st Shell 


2nd Shell 


ferns" ^ 






yr 


10 


0.21 


1.8 X 10^ 


,3.1 X 10^ 


100 


21 


1.8 X 10" 


■S.l X 10-* 


500 


105 


X 10-'' 


().2 X lO'"* 
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Fig. 1. — Three band false color image of the RCW 38 cluster using the IRAC bands 3.6 
yum in blue, 4.5 fim in green and 8.0 /um in red. The positions of the five identified bow 
shocks are indicated. 
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Fig. 2. — Enlarged view in three band false color IRAC/MIPS images of the five bow shocks 
in RCW 38. Above Left: The bow shock around YSO 293 in IRAC bands 3.6 fim, 4.5 /im 
and 8.0 fim. Above Right: The YSO 49 and bow shock in IRAC bands 3.6 /im, 4.5 fim and 
8.0 fim. Central Left: The YSO 129 in IRAC bands 4.5 /im and 8.0 /xm and MIPS 24/zm, 
with the bow shock only faintly visible in bands other than the 24 nm. The evaporating 
pillars are also bright at 8.0 and 24 /xm. Central Right: The fainter shock from YSO 803, 
to the SE in IRAC 4.5 /xm and 8.0 /xm and MIPS 24 /zm. It is only visible at 24/im. Below: 
The central shock from YSO 581 is shown in a three color image of IRAC 3.6 //m, 4.5 fim, 
8.0/im. IRSl glows brightly to the E of the bow shock, which is likely interacting with the 
strong winds from IRS2 further to the east in the image. 
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Fig. 3. — Spectral Energy Distributions of the five driving stars (stars symbols) and their 
shocks (square symbols). The 24//m fiuxes of YSOs 293 and 49 are contaminated by stellar 
emission. All fluxes for the YSO 581 shock are heavily contaminated by variable background 
and the relative location of 581 behind the shock. YSO 129 and source YSO 803 have only 
weak mid-IR fluxes and are identifled at 24//m. 
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Fig. 4. — Three band false colour image of RCW 38, with the smoothed ACIS-I 0.5-8 keV 
in blue, IRAC 4.5 fim in green and 8.0 /xm in red. Hot ionised gas outflows can be seen 
in the ACIS-I data. In the IRAC bands, ridges of emission from warm dust (8.0 fim) and 
shocked hydrogen (4.5 fim) coincide with the boundaries of the ionised gas outflows. Two 
of the stars with bow shocks, YSO 129 & YSO 293 (labelled a & b, respectively), also lie 
in regions of ionised gas, driven by winds from local massive stars, YSO 129 itself and two 
candidate O stars identified in Paper 1: ocl5-vl3 (c) and YSO 244 (d). 
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Fig. 5. — Grayscale of the difference between the 4.5 and 5.8//m IRAC image showing the 
central subcluster with ACIS-I data overlaid in contours. The subtracted image highlights 
the shocked hydrogen emission in the 4.5//m band in white. Three outflows can be observed, 
extending to the SW, NE, and NW in the ACIS-I contours. The class O/I (circles) and fiat 
spectrum (inverted triangles) protostars in the central core are overlaid showing the filament 
of possibly triggered star formation by the NW lobe of the outflow. Labels a-d as in Figure ID 
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Fig. 6. — Above: Grayscale of the 3.6/im IRAC image showing the location of DBS2003- 
124 with respect to the centre of ROW 38. Below: Three band false colour image of the 
newly identified NE subcluster in ROW 38, DBS2003-124, and the YSO-driven jet, with 
IRAC 3.6/im in blue, 5.8/im in green, and MIPS 24/im in red. The subcluster, in the upper 
left of the image, contains one source with bright 24/im emission, which is surrounded by 
tens of IRAC sources. Two candidate infrared HH-objects associated with the jet are visible 
in the centre and centre right of the image. The star powering the jet, NE21, lies NE of the 
brightest blue star in the image (OB candidate NE20). 



